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A power expansion scheme is set up to determine the Wigner function that satisfies the quantum 
kinetic equation for spin-1/2 charged fermions in a background electromagnetic field. Vector and 
axial-vector current induced by magnetic field and vorticity are obtained simultaneously from the 
Wigner function. The chiral magnetic and vortical effect and chiral anomaly are shown as natural 
consequences of the quantum kinetic equation. The axial-vector current induced by vorticity is 
argued to lead to a local polarization effect along the vorticity direction in heavy-ion collisions. 
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Introduction. — Chiral anomaly is an important quan- 
tum effect which is absent at the classical level. Recently 
it has been shown that such a microscopic quantum effect 
can have a macroscopic impact on the dynamics of rela- 
tivistic fluids, termed as the chiral magnetic and vortical 
effect (CME and CVE) [l-yl as manifested in currents 
induced by magnetic field and vorticity. Such effects and 
related topics have been investigated within a variety of 
approaches, such as AdS/CFT duality \^]^, relativistic 
hydrodynamics |9l4lll| , and quantum field theory [2|, Il2l - 
|l7|. However, it is still not clear how CME and CVE can 
emerge from a microscopic quantum kinetic theory. 

In this Letter we make a first attempt to derive both 
the CME and CVE from a quantum kinetic theory. A 
power expansion in space-time derivatives and weak ex- 
ternal fields is used to determine the analytic form of vec- 
tor and axial-vector components of the Wigner function 
that satisfies the quantum kinetic equation for spin-1/2 
massless fermions. The CME and CVE appear naturally 
in the induced currents. Chiral anomaly and other con- 
servation laws are also automatically satisfied. The axial- 
vector current induced by vorticity depends quadratically 
on the temperature, baryonic and chiral chemical poten- 
tial. So it should be present in both hot and dense mat- 
ter, and can lead to a local polarization effect in heavy-ion 
collisions as proposed in earlier studies [18-2(|. This pro- 
vides another possible future experimental measurement 
of the CVE in high-energy heavy-ion collisions. 

The quantum kinetic approach can provide a bridge 
between the microscopic and macroscopic description of 
the CME and CVE and should be more suitable for fu- 
ture simulations of both effects in heavy-ion collisions. 
The power expansion method can also be applied to the 
calculation of other transport coefficients. 



Quantum kinetic equation. — In a quantum kinetic 
theory, the classical phase-space distribution f{x,p) is 
replaced by the Wigner function W{x,p) in space-time x 
and 4-momentum p, defined as the ensemble average of 
the Wigner operator 2ll423J| for spin-1/2 fermions, 
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where "00, and '0/3 are Dirac spinor fields, x± 

are two space-time points centered at x with space-time 

separation y, and the gauge link U , 



U(x 
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ensures the gauge invariance of Wa/s . Here Q is the elec- 
tromagnetic charge of the fermions, and A^ is the electro- 
magnetic vector potential. Note that we use the metric 
convention g^" = diag(l, — 1, — 1, — 1). To simplify the 
quantum kinetic equation under a background field we 
consider a massless and coUisionless fermionic system in 
a constant external electromagnetic field F^j/ in the lab 
frame. Since we only consider a classical background 
field, we have dropped the path ordering in the gauge 
link in Eq. ^ . The Wigner function is a matrix in Dirac 
space and satisfies the quantum kinetic equation 2l|-|23l| , 



MP^ + T:^nw{x,p) = 0, 



(3) 



where 7'''s are Dirac matrices and V^ = d!^ — QP'^'-yd^. 
The Wigner function should contain information about 
quantum interactions and we will prove that all currents 
including chiral anomaly can be derived from the above 
equation. To this end, we decompose the Wigner func- 
tion in terms of 16 independent generators of the Clifford 



algebra, 



^{x,p) + ij'^:^{x,p) + 7^r^(a;,p) 



+7V-^M(a^.P)+2^'''^M-(^'P) 
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then leads to two decoupled sets of equations 
one of which relevant to our study reads, 



p^r^ = 0, p^i^^ 



0, 



v^-^, = 






0, V'^^^ = 0, 
-2{p^,s^^, ~p^£/^,), 



(5) 
(6) 

(7) 
(8) 



where e'"^'"^ is the Levi-Civita anti-symmetric tensor, 
'y^(x,p) and £/^{x,p) are the vector and axial- vector 
component of the Wigner function, which will give rise 
to the vector and axial- vector current, respectively, after 
integration over four-momentum. 

Power expansion. — We assume a system close to local 
equilibrium under a constant external field F*^'^ . There- 
fore, 'f'^{x,p) and s/^(x,p) will depend on x only through 
fluid four- velocity u(x), temperature T(x), chemical po- 
tential ijl{x) and chiral chemical potential iJi5(x). We will 
determine the analytic form of the Wigner function in 
terms of {p, F'^'^ , u, T, //, /us} from the kinetic equation. 

We further assume that the space-time derivative dx 
and the field strength F^v are small variables of the same 
order and can be used as parameters in the power ex- 
pansion of 1^ and J2^ (similar to the Knudsen number 
expansion in hydrodynamics). 
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where the subscripts 0, 1, ... denote orders of the power 
expansion. Note that 1^ and jz/^ are related to ^,f_i 
and 'y,^_i via Eqs. ([TllH]) [n > 1). One can therefore use 
an iterative scheme to solve 1^ and ^^ order by order. 

Note that the field strengths F'^'^ are assumed to be 
constant in the lab frame. Later, we have to define elec- 
tromagnetic fields in the local comoving frame of a fluid 
cell, E„ = uPF^p, B„ = (l/2)e„^,ypU^'F''P, which depend 
on X via the fluid velocity u{x). The space-time deriva- 
tive dx is then given by 



9f = d„T 
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Zeroth-order Wigner function. — In general, 1^'^ and 
£/q can only have two terms, each proportional to the 
zeroth-order four-vectors p'^ or u'^ with a total of four 
independent coefficients. Since the left-hand sides of 
Eqs. dlllHl) are at least of first order, the zeroth-order 
terms on the right-hand sides must vanish, which set the 
coefficients of the u^-terms to be zero. With additional 



constraints by Eq. ([5]), '^'^ and 
lowing forms. 



have to take the fol- 



y,^^pPS{p')Vo, <=p^<5(p^)Ao, (11) 

where Vq and Aq are the phase space distributions of 
massless spin- 1/2 fermions at the zeroth order and cannot 
be determined by Eqs. dMH]). We assume they take the 
equilibrium form, 

[Vo,Ao] - ^ 0(swp) [(/,,,? + /,,l), (/,,«- A,l)], 

s=±l 
f = ^ I 



,{X = R,L), (12) 



where R{L) denotes the right (left)-handed fermions and 
A*fl,L = fJ- ^ tJ-5 121 ■ Note that Vq (Aq) is the sum (dif- 
ference) of two positive distributions for any values of n 
and fj,^. This asymmetry between Vq and Aq as inputs to 
the iterative operation will feed down to the first-order 
Wigner functions Vf and Ai and the final vector and 
axial- vector currents, even though the kinetic equations 
in Eqs. ([SM]) are symmetric for V^ and A^. 

The zeroth-order Wigner functions should also satisfy 
Eq. (J6|), which provides constraints on fluid and thermo- 
dynamical variables. Substitute Eqs. (|imi2p into Eq. ([S|), 
we obtain V^'^'' and V^jz/q^ as sums of six indepen- 
dent terms involving the momentum vector p^ = ^apP'' 
(Ag-p = ga^p — UcrUp), tensor PaPp, scalars p^ and u ■ p. 
To ensure V^T^^ — V^jz/q'' = for any values of p, these 
six terms all have to vanish, resulting in the following 
constraints at the flrst order. 



A^^A"'' (daup + dpu^ - ^A^pAP^dpU,) 



TA'^Pdp!^ + QE'' = 0, 
u-du" - A'^PdplnT^Q, 



d^^^O, u'^d, ^ 



T 



T 



u-dT+ -TAP^dpU^ 



0, 
= 0. 



(13) 



Note that we have dropped S{po) terms from derivatives 
of 0{pq) and 0{—pq), which are irrelevant when carrying 
out the 4-momentum integration due to vanishing phase 
space at zero momentum. Since we are interested in cur- 
rents induced by external fields and vorticity, we consider 
only the static case with a constant temperature. The 
above constraints are reduced to. 



u-du" = 0, d„u" = 0, d„ii = -QE^ 
/i5 ~ const, (for T = const.) 



(14) 



which has a simple solution /i ~ const. — QE ■ x and a 
solenoidal fluid velocity u'^(x — ut/u^) with 9 • u = 0. 

First-order Wigner function. — With the zeroth-order 
Wigner functions in Eqs. (Illj|12p one can determine the 



first-order Y^^ and s/^ from Eqs. (I5][8]). A general form 
linear in the first-order variables X^ = (£''', B^,u!^) and 
constrained by Eq. ^ can be written as, 



lw^l 



J2 [uu{g''^-p''p^lp')p'{P ■ x)Zxi 

X=E,B,ui 



X3 



+e^^P''uxppX^Zx4] , 



(15) 



where ^^^ = {^^,£^1") and w^ = (l/2)e^^p<,u'^a''u'^ is 
the fluid vorticity. Note that X ■ u = 0. There arc 24 
independent coefficients Zxi — iVxi,Axi) [i = 1,2,3,4) 
in the above power expansion of the first order. With 
Eqs. pT|) and P^ for the zeroth and first-order Wigner 
functions, both of Eqs. (O and ^ at the first-order con- 
tain 3 different tensor structures, each consisting of terms 
linear in three independent variables X^^ = {E'^, B'^,uj'^). 
Setting these terms to vanish separately gives 18 equa- 
tions which leave only 6 of the 24 coefficients in Eq. ((T5]) 
undetermined. Further requiring Eq. ([6]) be satisfied by 
Yl^ and jz/f^, we can obtain the unique forms of 1^ and 
£/f_i to the first order, 

1 r)Z 

^^ = P^6{p^)Zo + -pAu''oj''-u-'^^]j^^6{p^) 

2 o(u ■ p) 



-Qp4u''B''~u''B'']ZoS'{p^) 
+Qe^'^P^uxppEaZoS'{p^), 



(16) 



where iT = (r, ^), Zq - (Fq, Aq) and Zo = (Aq, Vo). 

Induced currents, CME and CVE. — We can derive the 
vector and axial-vector current from the above Wigner 
functions up to the first order in power expansion. 



r 



■f-i- 
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= I d'^pY^' = nut'+£^Lot' +^BBt', (17) 

= f d^p£/^' ^ n^u^ + ^^u^ + iB5B'' ■ (18) 



The energy-momentum tensor T^'' can also be evaluated, 
T^" - i f d'^pip^'y+p'^r^') 

+ ^QC(^i''S" + u'^B^). (19) 

The charge n, n^ and energy density e in equilibrium. 

No = 27r J dpopl,[0ipo) - 0i-po)]ZNo, (20) 

are determined from the zeroth-order Wigner functions, 
where Nq = n^n^^e corresponding to i = 2,2,3 and 
Zno = Vo,Ao,Vb, respectively. The pressure is given 
by F = e/3. Coefficients ^, £_b, £.5 and ^55 are given by 

S = CTT / dpop'o[e{po) - e{-po)]Zso, (21) 



where 5 — £,,£b, £.5,^35 corresponding to j = 1,0,1,0, 
c = 2, Q, 2, Q, and Z^o = Aq, Aq, Vq, Vq, respectively. It 
is easy to verify following relations: ^ — {l/2)dnc,/dfj,, 
^5 = {l/2)dn/dfi, Cb - {Q/2)d£/dfi, and ^55 = 
{Q/2)dC5/dfM. 

One can complete the above integrals analytically to 
obtain coefficients £, £b, ^5 and £b5 of the induced cur- 
rents as functions of /i, /is and T , 
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i = ^^1^15, £b = ^/i5, 
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(22) 



^fi. (23) 



Thermodynamical quantities n, n^ and e can be similarly 
obtained. 

The current in Eq. (|17p induced by magnetic field and 
vorticity with coefficients (,b and ^ in Eq. ([2^ . known as 
the CME and CVE [2, 3, .9,], respectively, is a direct con- 
sequence of the quantum kinetic equation for the Wigner 
function. The axial-vector current in Eq. (IT51) induced 
by magnetic field and vorticity corresponds to some sort 
of reversed CME and CVE, respectively. These results 
are consistent to those obtained from the second law of 
thermodynamics in Refs. [10] and j24| except a quadratic 
term in temperature in ^5 induced by vorticity. It should 
be noted that Eqs. (J22ll23p . including the temperature 
term in ^5, have also been obtained independently in 
Ref. [17| within the Kubo formalism. 

Conservation equations for j^^ and j(^ 

d,r^O, d,j!^^^^E-B, (24) 

can be derived from Eqs. (J17m8p with constraints on fluid 
and thermodynamical variables in Eq. ([14]). The electric 
held in the chiral anomaly appears through dafJ- — —QE^ 
from Eq. (fT4| . Note that we derived the chiral anomaly 
here without regularization in contrast to the derivation 
in quantum held theory. This is because the Wigner func- 
tion contains two fermionic flelds separated in space-time 
(nonlocal) and therefore free of singularities. One can 
also verify the energy-momentum conservation equation 
in the background fleld. 



dpT'^'^ = QF-^^jp, 



(25) 



from Eqs. (ITTl) and p9| with constraints in Eq. (fT4|) . It is 
interesting to observe that constraints in Eq. (fT3|) or (fT4|) 
require uj^ \\ B^ \\ E^^, which is crucial for the energy- 
momentum conservation in Eq. ([25]). 

It is remarkable that we have derived from the quan- 
tum kinetic equation not only currents in Eqs. (|17m8p 
with their coefficients in Eqs. (|22ll23p but also a complete 
set of conservation equations with the chiral anomaly 
in Eqs. ()24ll25p for charge, chiral charge and energy- 
momentum, respectively. In contrast, these conservation 
equations are used as inputs to obtain the currents in 
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FIG. 1: (Color online) The axial current induced by vorticity 
leads to the local polarization effect. The momentum (spin) 
direction is in the red-dashed (blue-solid) arrow. 

Refs. [Sl llO| with the requirement of the second law of 
thermodynamics. 

Multi-flavor fluid. — So far we have only considered a 
fluid with a single type of fermions. An extension to the 
case of multi-flavor quarks is straightforward. We can 
consider a three-flavor fluid with u, d and s quark and 
their anti-quarks. Note that each quark carries Nc fun- 
damental color charges. For the induced electromagnetic 
and baryonic vector current j^ , 
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For this three- flavor quark matter we have ^^ Qf ~ 0, 
and ^^'^■■y"" = ^EM ^ Q^ rpj^jg implies that the CME 
(CVE) dominates the electromagnetic (baryonic) current 
[3|. For the induced baryonic axial- vector current jt^, 



^5 
C-B5 



i^'^S^f'^"^' 



Nc V- 



(27) 



Therefore, magnetic fields cannot induce the axial- vector 
current in a three-flavor quark matter, which can only be 
induced by vorticity. 

Local polarization effect. — An axial-vector current 
induced by vorticity implies that the right (left)-handed 
fermions move parallel (opposite) to the direction of vor- 
ticity. Since the momentum of a right (left)-handed mass- 
less fermion is parallel (opposite) to its spin, all spins are 
parallel to the direction of vorticity (see Fig. [T] for il- 
lustration). This results in the local polarization effect 
(LPE) similar to what was proposed in Refs. 1814201 due 
to spin-orbital coupling. The LPE can be measured via 
hadron (e.g. hyperon) polarization along the direction of 
vorticity or the global orbital angular momentum in non- 
central heavy- ion collisions [18]. Note that ^5 in Eq. ([?7)) 



has three quadratic terms in T, /x and /X5. Therefore, 
the LPE should be present in both high and low energy 
heavy-ion collisions with either low baryonic chemical po- 
tential and high temperature or vice versa. 
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